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We have studied the magnetism of the Pr^+ ions in PrFeAsO i- a; F^; [x = 0, 0.15) and its interaction 
with the Fe magnetic order (for x = Q). Specific heat data confirm the presence of a first excited 
crystal electric field (CEF) level around 3.5 meV in the undoped compound PrFeAsO. This finding 
is in agreement with recent neutron scattering experiments. The doped compound is found to have 
a much lower first CEF splitting of about 2.0 meV. The Pr ordering in PrFeAsO gives rise to large 
anomalies in the specific heat and the thermal expansion coefficient. In addition, a field-induced 
transition is found at low temperatures that is most pronounced for the magnetostriction coefficient. 
This transition, which is absent in the doped compound, is attributed to a reversal of the Fe spin 
canting as the antiferromagnetic Pr order is destroyed by the external magnetic field. 

PACS numbers: 75.40.Cx, 71.70.Ch, 75.80.+q, 65.40.De 



A. Introduction 

Layered FeAs-materials have been studied exten- 
sively since the discovery of superconductivity in 
LaFeAs0i_2;Fa; with transition temperatures Tc up to 
28 K. ^ By exchanging the nonmagnetic La with magnetic 
rare earths (RE) such as Pr or Sm the superconducting 
transition temperature could be increased above 50 K.-' "' 
All parent compounds of the i?£^FeAsO family {RE = 
La, Ce, Pr, Nd, Sm, Gd) behave rather similar They 
exhibit magnetic ground states with a structural distor- 
tion from a tetragonal to an orthorhombic lattice around 
150 K and subsequent formation of a spin density wave 
(SDW) of the Fe moments. Except for LaFeAsO addi- 
tional antiferromagnetic (AFM) ordering of the RE mo- 
ments is observed with transition temperatures of the 
order of 10 K. 

PrFeAsO has the highest RE ordering temperature 
among the i?£:FeAsO family with T^' w 12 K.'''^- The 
structural and SDW transitions take place around 153 K 
and 127 K, respectively. ' A delicate interaction between 
the Fc and Pr moments has been deduced from /iSR^'^ 
and Mossbauer ' ' experiments such that the ordering of 
each sublattice entails a reorientation of the other one. 
The magnetic properties of PrFeAsO are further compli- 
cated by crystal electric field (CEF) effects: In tetragonal 
symmetry the free-ion ground-state multiple! of Pr^+ is 
split into 5 singlets and 2 doublets, while the remain- 
ing degeneracy is lifted in any lower symmetry.'^ ' In- 
elastic neutron scattering (INS) experiments revealed a 
first CEF excitation around 3.5 meV that is split above 
T]^'. ' ' However, the energetic position of the other lev- 
els remained unclear. Substitution of O by F changes 
the magnetism fundamentally: It leads to a simultane- 
ous destruction of the Pr and Fe order. ' ' Moreover, addi- 
tional CEF excitations have been found around 10 meV 
for X = 0.13."' 



In this paper we investigate the Pr magnetism in 
PrFeAsO and its interplay with the Fe SDW. We start 
with an estimation of the CEF level scheme of Pr^"*" in 
PrFeAsOi_a;Fa; (x = 0, 0.15). Fluorine doping is found to 
have a significant influence on the splitting, in agreement 
with the INS results. We find an additional low-lying 
level around 2 mcV in the F-doped compound, which was 
not seen before. The sensitive dependence of the split- 
ting on the F substitution is discussed a possible rea- 
son for the absence of Pr magnetic order in the doped 
compound. Subsequently we present thermal expansion, 
magnetostriction and magnetization data. PrFcAsO is 
found to undergo a field- induced transition below T^^ , 
which is attributed to an Fe spin reorientation due to the 
destabilization of the AFM Pr order. 



B. Experimental details 

Polycrystalline samples have been prepared by solid 
state reaction as described in Ref. 18. The specific heat 
was studied in a Quantum Design PPMS by means of 
a relaxation technique. For the thermal expansion and 
magnetostriction measurements a capacitance dilatome- 
ter was utilized, which allows a very accurate study of 
sample length changes dL/L}'^ The linear thermal ex- 
pansion coefficient a was calculated as the first tempera- 
ture derivative of dL/ L, while the magnetostriction coef- 
ficient P is determined by the first derivative of dL/L 
with respect to the magnetic field B = iiqH. Mag- 
netization measurements were performed in a commer- 
cial VSM-Squid (Quantum Design) at an oscillation fre- 
quency of 14 Hz. 
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FIG. 1. (a) Specific heat Cp divided by temperature T for 
PrFeAsOi-xFi: {x — 0,0.15) in zero magnetic field and 9 T. 
For X = a clear anomaly with a maximum at 11.3 K is 
observed in zero field originating from the magnetic ordering 
of the Pr^"*" ions at T^"^. The specific heat of LaFeAsOo.gFo.i 
used as a phonon background is also shown for comparison, 
(b) Behavior of the anomaly at T^"^ of PrFeAsO in magnetic 
field, (c) The superconducting transition of the sample with 
X = 0.15 gives rise to a small anomaly seen in Acp/T = 
Cp(0)/T-Cp(9T)/T. 



C. The CEF splitting of Pr^+ 

First we will evaluate the crystal electric field (CEF) 
level scheme of the Pr'^"*" ions and the influence of fluorine 
doping by analyzing the speciflc heat of both compounds. 

Fig. la shows the speciflc heat Cp of PrFeAsO i_a;Fa; 
{x = 0,0.15) as Cp/T for T < 60 K. Data for a; = at 
higher T are similar to those published in Ref. 20 with 
anomalies due to the structural transition at Tg — 145 K 
and the SDW formation at = 129 K.-^ These anoma- 
lies are absent for x = 0.15. Here, we focus on the 
regime below 60 K. The undoped sample exhibits a sharp 
anomaly with a maximum at 11.3 K, which is on top of a 
broad hump. This behavior is attributed to the magnetic 
ordering of the Pr"^"*" moments in the presence of CEF 
splitting. The evolution of the anomaly in an applied 



magnetic field is shown in Fig. lb. Small fields of 3 T 
and 5 T lead to a weak shift to lower T and a reduction 
in height. In 7 T, the anomaly is hardly seen anymore. 
Application of a field of 9 T suppresses the anomaly com- 
pletely and leaves only the hump and a tiny kink. This 
behavior is in line with a suppression of the antiferro- 
magnetic Pr order by the external magnetic field. On 
the other hand, the hump due to the CEF splitting is 
scarcely changed by the field, because the corresponding 
energy scale is much larger than the Zeeman splitting. 
For the sample with 15 % fiuorine doping, no Pr order- 
ing is observed down to 1.5 K, and the hump is shifted 
to lower T. In addition, a small anomaly is found at the 
superconducting transition, which is largely suppressed 
in 9 T. The difference between the zero-field and 9 T 
data reveals a small A-shaped anomaly at Tc = 42 K 
(cf. Fig. Ic). 

There is a significant magnetic contribution c,nag to the 
specific heat of PrFeAsOi-^^Fa, stemming from the Pr"^"*" 
ions. An exact determination of Cmag is impossible due to 
the presence of the various phase transitions. In order to 
get a rough estimate of the magnetic contributions c^^g, 
we used LaFeAsOo.gFo.i as a nonmagnetic reference. 
This is reasonable as substitution of the rare-earth ion is 
not expected to change the phonon spectrum significantly 
and the anomaly due to the superconducting transition of 
LaFeAsOo.gFo.i at 26 K is very small. Therefore, its con- 
tribution to Cp is negligible compared to the Pr contribu- 
tion. However, one has to keep in mind, that there might 
be a different electronic contribution to Cp (see discus- 
sion below). cj^*g is obtained by subtracting the specific 
heat of LaFeAsOo.gFo.i, which is shown for comparison 
in Fig. la, from the one of PrFeAsOi-^jFj;. The result is 
plotted in Fig. 2 as c^ljT. Both PrFeAsOi-^^F^j, sam- 
ples exhibit Schottky-like anomalies, whereas the maxi- 
mum occurs at lower T for the sample with fluorine dop- 
ing. This demonstrates already that the doping changes 
the CEF level scheme and results in a lowering of the flrst 
excited level. 

For PrFeAsO0.85F0.15 no magnetic ordering is observed 
down to 1.5 K and Cmag is only due to the thermal oc- 
cupation of higher CEF states. We may use cj^^^g to es- 
timate the CEF splitting. The data range between 20 K 
and 45 K, which contains the superconducting anoma- 
lies, was omitted for the fitting. The best description of 
the data is obtained for a ground-state doublet (D) with 
excited singlets (S) at 26 K and 40 K. Two more lev- 
els are found at 140 K (quasi-triplet (T)) and at 400 K 
(quasi-doublet). These two levels may in fact consist of 
three or two close-lying singlets, which, however, cannot 
be resolved from our data. In our fit we also accounted 
for a difference between the electronic contributions of 
LaFeAsOo.gFo.i and PrFeAsO0.85F0.15. Below the super- 
conducting transitions it was taken zero assuming that 
all charge carriers are condensed to Cooper pairs. Above 
Tc wc obtain a value of A7 = 25 mJ/(molK^) in ad- 
dition to the Sommerfcld coefficient of LaFeAsOo.gFo.i 
of 3.5 mJ/(molK^).-' Thus, we may estimate 7 = 
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FIG. 2. Estimate of the magnetic contribution to the specific 
heat of PrFeAsOi-iFa; {x — 0, 0.15) calculated by subtracting 
the data of LaFeAsOo.gFo.i. The sohd hues are the results 
of a fit to a Schottky contribution from CEF spUtting for 
X — 0.15. Above Tc an additional charge carrier contribution 
A7 = 25.2 mj/mol K is taken into account, which leads to 
a shift of the line. The dotted hne is the calculated Schot- 
tky contribution for a ground state singlet and an excited 
doublet at 41 K. The agreement with the data for a: = is 
much improved by adding a second doublet at 160 K (dashed 
line). The inset shows the estimated contribution to Cmag 
of a; = 0.15 from sites with one or more fluorine neighbors 
in a model with different CEF schemes for differing environ- 
ments^'' as explained in the text. The dotted line is the cal- 
culated contribution expected for two singlet states at and 
23 K for 37% of the Pr sites. 



31.5 mJ/(molK^) for PrFeAsO0.85F0.15, which is com- 
parable to the value of 34 mJ/(molK^) determined for 
PrFeAsO.20 

The result of the fitting is shown as lines in Fig. 2 
together with a schematic drawing of the correspond- 
ing level scheme. Our fit describes the data very 
well. However, the proposed ground-state doublet raises 
the question, why no magnetic ordering is observed in 
PrFeAsO0.85F0.15- We would like to mention, that a de- 
scription of the data using a ground-state singlet was 
not possible. This might indicate that c^*g is not sufB- 
ciently precise to allow estimating the CEF contribution. 
However, the presence of two singlets at 26 K and 40 K 
above a ground-state doublet suggested by our fit is de- 
termined mainly by the data below 20 K. In this temper- 
ature range, the estimate for Cmag is rather good, because 
the phonon contribution is comparably small. A possible 
answer is given by a recent neutron diffraction experi- 
ment. Goremychkin et al. studied CEF excitations in 
PrFeAsOi_^F^ (x = 0,0.13)."' The undoped compound 
was found to have an excitation around 3.5 meV that 
is split above Tj^' . Substitution of fluorine removed the 
splitting and led to a reduction of that peak. In addition, 
two more peaks were observed at 9.7 meV and 11.8 meV. 



The authors explained their observation by the presence 
of two different well-defined charge environments result- 
ing from a random distribution of fluorine on the oxygen 
sites. This leads to 5 different nearest-neighbor configu- 
rations for the Pr'^+ ions, the most common of which are 
those with none (58 %) and one (35 %) F neighbor. The 
excitation at 3.5 meV with reduced height was attributed 
to the Pr'^^ ions with unchanged nearest-neighbor config- 
uration, while the two additional peaks were attributed 
to Pr"^"*" ions with one fluorine as nearest neighbor. 

This idea can be also applied to our specific heat data: 
We assume that cj^^g of Pi'FeAsOo.85Fo.i5 contains a 
CEF contribution from the 52 % Pr^+ ions (for x = 0.15) 
without fluorine neighbor. It can be estimated from c°^*g 
of PrFeAsO, because the entropy change related to the 
magnetic ordering in the undoped compound is rather 
small. After subtracting this contribution we end with 
a magnetic specific heat Cmag-iF dominated by the CEF 
contribution from the 37 % ions with one F neighbor. 
The resulting curve is shown in the inset of Fig. 2. The 
low-temperature part is dominated by a Schottky-like 
anomaly, which is surprisingly well described by two sin- 
glets at and 23 K at 37 % of the Pr sites. Due to 
the various approximations used so far, we refrain from 
a further analysis of Cmag-iF- However, the broad max- 
imum observed at higher temperatures may be ascribed 
to the CEF levels around 10 meV found in the neutron 
scattering experiments. The splitting of 23 K determined 
from Cmag- IF is close to the lowest splitting of 26 K de- 
termined by our fit of c^ag- Therefore, our data clearly 
prove the presence of a first excited level at about 2 meV 
for PrFeAsO0.85F0.15- No respective excitation was ob- 
served in the neutron scattering experiment. We suggest, 
that it probably corresponds to a forbidden transition. 

We now turn to the undoped compound. A CEF ex- 
citation around 3.5 meV has been found in the neutron 
scattering experiment mentioned above. We use this 
value to model our magnetic specific heat. A quite good 
description of the low-temperature part is indeed ob- 
tained assuming a singlet ground-state with a doublet 
at 41 K corresponding to 3-5 meV (cf. dotted green line 
in Fig. 2). The agreement at high T is significantly im- 
proved by assuming a second doublet at 160 K, which is 
at the limit of the measurement range in Ref. IG- The 
other states are supposed to lie at even higher energy. 

Summarizing the analysis presented so far we conclude 
the following: Our specific heat data for PrFeAsO are 
consistent with a singlet ground state and excited dou- 
blets around 41 K and 160 K, whereas the remaining 
levels lie at higher energy. PrFeAsO0.85F0.15 has a first 
excited state at significantly lower energy corresponding 
to about 23-26 K. This low-lying singlet is responsible 
for the shift of the hump in Cp to lower T (cf. Fig la)- 
Our data are also in line with neutron scattering results, 
which suggest different Pr sites for PrFeAsO0.85F0.15, de- 
pending on their environment. In this model, 52 % of the 
Pr ions have no fiuorine neighbor and CEF levels similar 
to PrFeAsO with a first excited state being a doublet at 
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FIG. 3. Temperature dependence of the linear thermal ex- 
pansion coefficient a(T) of PrFeAsOi-iF^ {x = 0,0.15) in 
magnetic fields of T and 9 T. Clear anomalies are seen in 
the zero-field thermal expansion curve for a; = at the struc- 
tural/SDW transitions and at the magnetic ordering tem- 



perature of the Pr 



These are absent for the doped 



compound. Application of a magnetic field leads to a strong 
change of a{T) for PrFeAsO below 70 K. By contrast, there is 
almost no effect on a(T) for the doped compound. The inset 
shows the low-T part on a larger scale. 



41 K. The 37 % Pr ions with one fluorine neighbor have 
a different splitting, whereas the first excited state is a 
singlet at 23 K. 

The change of the CEF level scheme upon F substitu- 
tion is probably also responsible for the absence of the 
Pr magnetic order in PrFeAsOo.ssFo.is- Another possible 
reason might be the absence of the Fe SDW. In fact, it has 
been found that both F doping and Ru substitution on 
the Fe site lead to a concomitant suppression of T^'^ and 
jiPr I . ,22 -pjj^jg suggests that both orderings are linked. 
However, competing magnetic structures have been pro- 
posed for the Fe and Pr sublattices, ' ' as explained in 
more detail below. This renders such a connection rather 
unlikely. Instead, in view of the sensitive dependence of 
the first excited CEF level on the Pr environment we sug- 
gest that Pr ordering takes place only for a specific CEF 
splitting. One may even speculate that the "CEF dis- 
order" due to the presence of different Pr environments 
and consequently different CEF states is responsible for 
the absence of Pr magnetic order in PrFeAsO0.85F0.15- 



D. The interplay of Fe and Pr magnetism 

We now turn to the interplay of Fe and Pr magnetism 
in PrFeAsO. For this purpose we present thermal expan- 
sion, magnetostriction, and magnetization data. 

The linear thermal expansion coefficient a of 
PrFcAsOi-ajF^; {x = 0,0.15) measured in zero magnetic 



field and 9 T is shown in Fig. 3. The structural and 
SDW transitions of the undoped compound give rise to 
a large, broad anomaly around 150 K.^'^ Both transitions 
cannot be distinguished in our data, probably due to the 
small height of one of the anomalies. The ordering of the 
Pj.3+ jnomcnts in PrFcAsO shows up as another, posi- 
tive anomaly with a maximum at 11.4 K, close to the 
maximum in Cp/T. This anomaly is rather sharp and 
symmetric, which suggests that the transition might in 
fact be of first-order type. The shape of our specific-heat 
anomaly is also in line with a broadened first-order tran- 
sition, yet we cannot exclude, that it is a second-order 
one. 

Although the pronounced background hinders the pre- 
cise determination of the anomalous changes, the hydro- 
static pressure dependence can be extracted quantita- 
tively from the anomalies in a and Cp by means of ei- 
ther the Clausius-Clapeyron or the Ehrenfest relation 
depending on the nature of the phase transition. Sup- 
posing a weak first-order character of the anomaly, ana- 
lyzing the data in Fig. lb and 3 yields a volume change 
of Ay / 3 ■ Aa{T)dT w 3.9 ■ lO'^/K and an entropy 
change of AS" w 0.13 J/(mol K) at the transition. Ap- 
plying these estimates we obtain the hydrostatic pressure 
dependence: 
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1.2K/GPa. 



(1) 



A similar value is obtained assuming that the transition 
is of second order. In this case the pressure dependence 
is determined from the jump heights Aa 1.0 • 10~^/K 
and Acp w 1.4 J/(mol K), which yields: 



"^-^ N 

dp 



TV, 



3Aa 



mol 



l.OK/GPa. 



(2) 



So far, no pressure experiments have been performed on 
PrFeAsO. However, the estimated pressure dependence 
of r^' is comparable to the one of the Ce ordering tem- 
perature in CeFeAsO of 0.9 GPa/K.- ' 

A field of 9 T leads to a strong decrease of the thermal 
expansion coefficient of PrFeAsO below approximately 
50 K. The sharp anomaly at 11 K is suppressed. Instead 
a broad maximum is observed around 8 K as seen in the 
inset of Fig. 3, which shows the low-T part of a(T) on 
a larger scale. Since no anomaly is found in the corre- 
sponding specific heat curve, this feature is not related to 
a phase transition, but rather to thermal population of 
higher Pr'^"'" states. Contrary to the undoped compound, 
the thermal expansion coefficient of PrFcAsO0.85F0.15 ex- 
hibits a smooth temperature dependence. No anomaly is 
found at the superconducting transition similar to the 
findings in LaFcAsOi-^jFa;. ' ' Most probably, it is too 
small to be seen in our data. Application of a magnetic 
field of 9 T leads to a lowering of a, which, however, is 
much weaker than for the undoped system. 

In view of the large change of a of PrFeAsO induced 
by a magnetic field we expect a sizeable magnetostrictive 



5 



2.0 



"? 1.5 
o 



T3 



II 



1.0 



0.5 



0.0 



PrFeAsO, F 

1-X X 




X = 0, r= 5 K>, 






y ^^"^^ ^"^^15K" 


7.5 


r ^^^K" "20 K. 


10 K 


\-^30 K" 








x = 0.15,T=5K \ 

.linn k- 





5 10 

6(T) 



15 




e(T) 



FIG. 4. Magnetostriction coefficient /? vs. magnetic field B 
of PrFeAsO at different temperatures up to f 00 K (full sym- 
bols). For T = 5 K, 7.5 K, and 10 K kinks marked by arrows 
are observed at 5.8 T, 5.4 T, and 3.7 T, respectively. At 
higfier T the transition disappears. For comparison data for 
PrFeAsO0.85F0.15 measured at 5 K are shown (open symbols). 
At higher temperature, /3 is even smaller for this compound. 



FIG. 5. The inset shows an example for the field dependence 
of the magnetization M{B) of PrFeAsO at 10 K. M{B) ex- 
hibits an almost linear behavior with small deviations above 
approximately 4 T, as illustrated by the dashed line. They 
are better seen as a step in dM/d_B shown in the main plot 
for different temperatures. The solid lines are a fits to an 
empiric formula as explained in the text. For clarity, not all 
investigated temperatures are shown. 



effect in the material. In Fig. 4 we present the mag- 
netostriction coefficient /? = 1/LdL/dB of PrFcAsO at 
different temperatures up to 100 K and in magnetic fields 
up to 16 T. At high temperatures, an almost linear de- 
pendence l3{B) is found. With decreasing T, the magne- 
tostriction coefficient increases, and a maximum in I3{B) 
develops at higher fields. Below the Pr ordering tem- 
perature, the behavior of (3{B) changes qualitatively. In 
low fields, 13(B) increases linearly with a temperature- 
independent slope. Interestingly, this slope is smaller, 
than the one observed at 15 K, i.e. above the phase tran- 
sition. At higher fields, kinks are observed in [3{B) at 
5.8 T, 5.4 T, and 3.7 T for 5 K, 7.5 K and 10 K, respec- 
tively. For fields larger than approximately 10 T, the 
low-T magnetostriction resembles the field dependence 
observed above 15 K. For comparison, a curve measured 
at 5 K on PrFcAsOo.s.sFo.is is also shown in Fig. 4. As 
expected from the small change of the thermal expan- 
sion coefficient in magnetic field, the magnetostriction is 
much weaker than for the undoped compound. 

The magnetostrictive effect observed in PrFeAsO be- 
low about 50 K, i. e. well above T^" is rather large. 
This is surprising since the magnetostrictive effect for 
PrFcAsO0.85F0.15 is significantly weaker. Apart from the 
Pr magnetic ordering there are two important differences 
between the undoped and the doped compound: (1) In 
addition to the Pr order, the Fe spins form a SDW below 
129 K in PrFeAsO. It may not be directly responsible for 
the large magnetostriction. For comparison: measure- 
ments on LaFeAsO allow estimating an upper limit of 
/3 < 5 X 10~*/T for this compound, despite the presence 
of Fe SDW order.-" However, the Fe SDW induces a po- 



larization of the Pr'^^ moments already well above Tj^' .^'^ 
The applied magnetic field thus acts on a compound with 
net magnetic moments from both the Fe SDW and the Pr 
sublattice, which may be the reason for the large magne- 
tostrictive effect. (2) As evident from neutron scattering 
data"' and our specific heat analysis, fluorine substitu- 
tion changes the CEF level scheme of part of the ions. 
This may explain the observed suppression of l3{B) at 
least to some extent. Only 52 % Pr ions without fluo- 
rine neighbor have an unchanged CEF scheme. There- 
fore, one expects about half the magnetostrictive effect 
of PrFeAsO for the fluorine-doped compound from these 
sites. Moreover, wc do not know the magnetostriction 
contribution of the remaining 48 % Pr ions with at least 
one fluorine neighbor. In addition, the strain induced by 
the differing favored expansion coefficients of neighboring 
cells will lead to a more difficult behavior than a simple 
sum of the effects. Measurements on samples with dif- 
ferent fluorine content might clarify the relevance of the 
specific CEF level scheme and the presence of the Fe 
SDW for the large magnetostriction of PrFeAsO. 

Now wc turn to the transition observed in the magne- 
tostriction coefficient of PrFeAsO. It is expected to give 
also rise to an anomaly in the magnetization M{B) of 
the material. We performed measurements of M{B) in 
fields up to 7 T and for temperatures between 1.8 K 
and 50 K. The magnetization is found to increase al- 
most linearly with field. A small change in slope is found 
for T < 12.5 K, that is more obvious from the deriva- 
tive AM /c\B. As an example, the inset of Fig. 5 shows 
M(B) measured in 10 K. Deviations from a linear be- 
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FIG. 6. The magnetic phase diagram of PrFeAsO, which was 
obtained from the magnetostriction (green triangles) and the 
magnetization (black squares) data. The error bars for the 
magnetization results have been determined from the width 
of the transition in dM/dB as explained in the text. The red 
crosses indicate the positions of the maxima in Cp/T, while 
the blue triangle corresponds to the maximum in the thermal 
expansion coefficient. The adjacent dotted lines mark the po- 
tential range for the phase transitions. The zero-field ranges 
for Cp/T and a are almost identical and therefore hard to dis- 
tinguish in the plot. The dashed line is a guide to the eye. 



havior arc observed above 4 T as indicated by the dashed 
line. This corresponds to a step at this field in the deriva- 
tive dM I AB shown in the main plot of Fig. 5 for different 
temperatures. With increasing temperature, the transi- 
tion field i?o shifts to lower B. At 12.5 K the transition 
has disappeared. In order to determine the position of 
the transition we fitted AMjdB with an empiric formula 
dM/dS = Ao AiB - hl\\ ^ [B/BqY]. For z > 2 this 
equation describes a broadened step-like function with 
step hight h. The broadening is determined by the pa- 
rameter z. At Bq half of the step height is reached. The 
parameters Aq and Ai account for a linear background. 
As a measure for the uncertainty in Bq we take the full 
width at half maximum (FWHM) of the derivative of the 
fit. The fits are also shown in the main plot of Fig. 5. 

Taking the transition fields determined from the mag- 
netostriction coefficient and the Bq values determined 
from the magnetization curves, we can draw a magnetic 
phase diagram for PrFeAsO, which is shown in Fig. 6. 
The error bars given for Bq are determined from the 
FWHM in the derivatives of the fits as explained above. 
Unfortunately, our specific heat and thermal expansion 
data do not allow extracting clear transition tempera- 
tures. In particular, no unambiguous conclusion on the 
character of the transition even in zero and small mag- 
netic fields is possible since it can be either continuous 
or weakly first-order. However, the transition is clearly 



suppressed in external magnetic field and the peak max- 
ima provide a lower limit for the transition temperatures. 
The actual transitions may occur between the maxima in 
Cp/T{T) and a(T) and the kinks in the data at the end of 
the anomalies. These ranges are marked by dotted lines 
in Fig. 6. In addition we also plot the peak maxima of 
the anomalies. 

Fig. 6 evidences what has been noted before: The kinks 
in M{B) and P{B) are observed only below T^''. In fact, 
the maxima of the anomalies in Cp/T{T) and a{T) in 
zero field fit well to the phase line. However, the specific- 
heat anomaly in 5 T and most probably also the one in 
3 T is observed at higher T. More precisely: the line 
marking the possible range for the Pr magnetic ordering 
transition in 5 T clearly does not cut the phase border. 
This suggests, that the transition observed in M{B) and 
f3{B) has a different origin and is not caused directly by 
the field-induced suppression of the Pr order. In fact, 
the observation of a kink in both properties is a rather 
unusual behavior for such an effect. Typically an abrupt 
reorientation of magnetic moments shows up as step-like 
anomalies in 13(B) and M{B). 

One may speculate that the presence of two magnetic 
subsystems with competing magnetic structures is re- 
sponsible for the observed behavior. A model for the 
magnetic structure of PrFeAsO has been proposed from 
/iSR experiments:^ ' Below 127 K, the Fe moments order 
in an antiferromagnetic stripe structure parallel a:' The 
resulting internal field induces a magnetic moment on the 
Pr'^+ site along c. The Pr order below Tj^'' leads to a re- 
orientation of both subsystems. The Pr"^"*" moments are 
now oriented along a, which, in turn, leads to a canting 
of the iron moments with a significant component along 
c. A Fe spin reorientation upon Pr ordering was also de- 
duced from Mossbauer spectroscopy.^^ Thus, it appears, 
that the preferred magnetic structures of the Pr and Fe 
subsystems compete with each other. 

We suggest the following picture for PrFeAsO: In zero 
field the Pr moments order below T^'^, which is accom- 
panied by a canting of the Fe spins. Application of a 
magnetic field destabilizes the antiferromagnetic Pr or- 
der as confirmed by the observed small shift of the spe- 
cific heat anomaly to lower T. This weakening of the Pr 
order allows the Fe sublattice to keep its preferred ori- 
entation down to lower T. Above a field of about 6 T 
the Pr ordering is suppressed completely and only a po- 
larization by the Fe SDW and the external field remains. 
With regard to our field-dependent measurements, the 
weakening of the Pr order leads to a decoupling of the 
Fc SDW above a certain field, at which the Fe spins re- 
sume their preferred orientation. The kink observed in 
our magnetostriction data then does not correspond to 
a sudden flip of Pr moments upon increasing field, but 
rather to a slow rotation, which changes speed, as the Fe 
SDW is released. This picture of a field-induced rotation 
of the Pr moments accompanied by a decoupling of the 
Fe SDW is also supported by recent measurements of the 
resistivity p."" A weak maximum was observed in p(T) 
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around 6 K in fields up to 6 T. Since the the FeAs layers 
are responsible for the charge transport it was attributed 
to a reorientation of the Fe moments induced by the Pr 
order below T^''. The authors corroborated their inter- 
pretation by the observation of a broad hump around 
6 K in the derivative of the specific heat divided by tem- 
perature. This latter feature may be also ascribed to 
the Schottky anomaly from the excited doublet at 41 K, 
which has a maximum in (l{c-axa.glT) / dT at 6.7 K. On 
the other hand, the interpretation of the resistivity max- 
imum as a signal of Fe moment reorientation caused by 
the Pr order is in line with our discussion: The maxi- 
mum in p{T) is observed only in fields below 6 T. At 
higher fields the Pr order is very weak or even destroyed 
as confirmed by the tiny anomaly in Cp observed in 7 T 
(sec Fig. lb) . As a result the Fe moments keep their pre- 
ferred orientation and the maximum in p{T) disappears. 
Moreover, a kink was found in p{B) at low temperatures 
by several groups.""'"'''^' The positions of these kinks fit 
rather well to our phase line, whereas small deviations 
may be due to the slightly different values. There- 
fore, we suggest that these kinks are also a signature of 
the field-induced rotation of the Pr ions accompanied by 
Fc spin reorientation. 



E. Summary 

In summary we have studied the Pr magnetism in 
PrFcAsOi_;j;F^, [x = 0,0.15). The CEF level scheme of 
the Pr'^"'" ions is strongly infiuenced by fluorine doping. 
The first excited CEF level in PrFeAsO is found around 
40 K, whereas a significantly lower splitting correspond- 
ing to about 23-26 K is determined for PrFeAsOo.ssFo.is- 
This lowering as well as a different CEF splitting for Pr 
sites with different environment may be the reason for 
the absence of Pr order in the fluorine-doped compound. 

At low temperatures, a fleld-induced transition is 
found for PrFeAsO, that is not directly related to the 
suppression of the AFM Pr order. Instead it comes from 
the interplay between the Fe and Pr moments, which 
appears to be very sensitive to application of magnetic 
fields. The transition is attributed to a reversal of the 
Fe moment canting induced by the Pr ordering as the Pr 
order is suppressed in magnetic field. 
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